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Abstract 



By designating one eigenvector of the mass matrix, one can reduce the free parameters in 
the mass matrix effectively. Applying this method to the quark mass matrix and to the lepton 
mass matrix, we find that this method is consistent with available experimental data. This 
approach may provide some hints for constructing theoretical models. Especially, in the lepton 
^is^ . sector, the Koide's mass relation is connected to the element of the tribimaximal matrix through 

' Foot's geometrical interpretation. In the quark sector, we suggest another mass formula and 

, the same procedure also applies. 

^ . I. INTRODUCTION 

^ i' The understanding of quark masses and mixings has posed a major challenge in particle physics 

for a long time. Recently, the non-zero neutrino masses and their mixings have been confirmed [T], 
^ I which implies that the mixing also exists in the lepton sector, just like that in the quark sector. 

<^ . A key step to understand the masses and mixings of quarks and leptons is to determine the mass 

matrices of quarks and leptons. One popular method, suggested by Fritzsch [5], is the texture 
zero structure. For example, the four texture zero structure can survive current experimental 
' tests [3]. In the lepton sector, other matrices, for example, based on the Vfj.-Vr symmetry and the 

A4 symmetry, have been suggested [T]. Especially, the nearest- neighbor- interaction form (NNI- 
form), can be implemented in some grand unified theories [4], and is consistent with the current 
experimental data from quarks and leptons [U |3] . Although the progress have been made in these 
directions, there is still no a commonly granted standard theoretical model for these problems. 
Therefore, other phenomenological approaches are necessary and worthy to be explored, which may 
provide some hints for constructing theoretical models. In this paper, we explore a way that can 
realize Koide's mass formula [6] through Foot's geometrical interpretation [7]. 

The Yukawa sector of the standard model has too many free parameters. In order to make 
definite predictions, we must make efforts to reduce the redundant parameters effectively. As we 
have emphasized, many papers have been devoted for this purpose. One common character of 
these papers is to reduce the redundant parameters by virtue of some symmetry [U H]. In this 
paper, we explore another way, which is different from these approaches. The main ideas are as 
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follows. First, in the standard model, we can choose the mass matrices to be Hermitian [81. Then by 
designating one eigenvector of a mass matrix, we can reduce the redundant parameters effectively. 
In fact, as we will illuminated below, if we make some assumptions and input the values of the 
mass parameters, only four free parameters are left. It is well known that the Cabibbo-Kobayashi- 
Maskawa (CKM) [9] matrix has four parameters. Therefore, adjusting the values of these left free 
parameters, we can fit the experimental data in principle. We find that this way is consistent 
with available experimental data. Because in this approach we have the freedom to choose the 
eigenvectors, it gives us some advantage to realize extra goals. For example, in the lepton sector, 
Koide's mass formula is connected to the entry 33 of Maki-Nakagawa-Sakata (MNS) matrix [TU] 
through Foot's geometrical interpretation, and it is consistent with available experimental data. In 
the quark sector, we suggested another mass formula, which can be connected to the entry 33 of 
CKM matrix through Foot's geometrical interpretation. 

This paper is composed of five sections. In Sec. II, we introduce the method in detail. In Sec. Ill 
and IV, we apply this method to the quark sector and to the lepton sector respectively. We make 
some conclusions in Sec.V. 



II. THE METHOD 

In the standard model, the mass matrices are complex matrices in general, but we can use the 
freedom of right-hand rotation to make them Hermitian 8J . So without loss of generality, we start 
our discussion from Hcrmite mass matrices. 

Supposed a general Hermite matrix 
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in which A, B, C , (pD, (f'E, (f'F are real and Z?, E, F are nonnegative. 
The matrix M can be written in another way 

_ A F D 

M = P^AfP = ptj p B Eexp'" \ P, (2) 
D E exp-*" C 

in which P — diag(l, exp*"^*" , exp"^° ), and a = (pE — <j)D + 4>f- Given that M has one eigenvector 
~x'^ = (xi, X2 exp'*^, ^3 exp*''')"'" belonging to its eigenvalue A, we have the eigenequation 

Ulc = A^, (3) 

in which xi, X2 and are nonnegative real numbers; (3 and 7 are real numbers. 
After some simplification, it reduces to 

Xi 

a;2exp*(^+'^^^ ) =0. (4) 
a;3 exp*(T+"*°) 

We see that Eq. (4) contains complex variables, and it will be difficult to solve them. We notice 
that it will be simple in a special case, in which we let /3 = — 0f and 7 = —4>d- This implies that we 
designate a real eigenvector to the matrix M in Eq. By this choice, all of the equations have real 
variables, and they can be solved with little labor. It is obviously that this is only a conventional 
choice, which simplifies the equation effectively. Of course, we should consider the possibility that 
this choice may be not appropriate, hence the equations have no solutions. However, in Sec. Ill 
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and Sec. IV, we will give the numerical results, which imply that our choice is compatible with 
experimental data. In this paper, we will restrict our discussions on this simple case. Of course, 
other cases, in which (3 + (ftp 7^ and 7 + (pD 7^ 0, are not excluded if they are needed. Then we 
have 

( A-\ F D \ ( xx\ ( xx\ 

- 0, (5) 

in which / is the identity matrix. In Eq. ([5]), as A, B, C, D, F, A, xi, xi and x-j, are real 
numbers, if E and 3:3 are nonzero, a must equals to or tt. Therefore, this matrix identity produces 
three equations. It is well known that two of these equations are independent with each other. We 
choose the independent equations to be 

{A~\)xx+Fx2 + Dxt. = 0, (6) 
Fxi + (B-A)a;2 + £:exp^"x3 = 0, (7) 

in which a = or tt. In addition to Eq. JH) and Eq. d?]), we have three eigenequations 

{A - Xi){B - Ai)(C - Ai) - E^(A - Ai) - D^(B - Ai) - F'^{C - Ai) + 2DEFN = 0, (8) 
(A - X2)iB - X2){C - X2) - E^iA - Ai) -D\B- X2) - F\C - X2) + 2DEFN = 0, (9) 
{A - A3)(B - Xi){C - A3) - E'^{A - A3) - - A3) - - A3) + 2DEFN = 0, (10) 

in which N — cos a = ±1. 

So far, we have five equations, and we have six free parameters, among which only one parameter 
is still free. We can let F to be the free parameter. Once we fix the value of F, all other parameters 
are fixed. These equations can be solved analytically, but the expressions are too complicated. In 
order to simplify the expressions of the solutions, we give some analysis in Appendix B. When 
we apply them to the lepton sector in Sec. Ill and to the quark sector in Sec. IV, we will give the 
analytical expressions explicitly in Appendix C and in Appendix D. However, when we adjust the 
left free parameters to fit the experimental data, the numerical approach is needed. In the text, 
we display the numerical results. Also, it is possible that these equations have no solutions, if the 
eigenvalue and the eigenvector are not appropriate. However, in Sec. Ill and Sec. IV, we will show 
that for the parameters we choose, the solutions always exist, as we will display explicitly. 

The key point of our method is to choose the appropriate eigenvalue and the appropriate eigen- 
vector. In the following application, we will choose the eigenvector and eigenvalue according to 
physical ground. 

With the method we suggested above, if we fix the value of the left free parameter, we can 
fix the matrix M . Now we turn to show how we can use this method to determine the mixing 
matrix. We take the CKM matrix for example. We let the up-quarks mass matrix to be M„, and 
the down-quarks mass matrix to be Md- Like M, we can write M„ and Md as 

= PiM^Pu, M~d = P^MdPd. (11) 

We designate 'y''^ = (yi, 2/2, J/a)^ as the eigenvector of M„ belonging to its eigenvalue A^, and 
'^'^ = {zi, Z2, 2:3)^ as the eigenvector of Md belonging to its eigenvalue A^. Then we have 

(M„ - A„/) I L 1 = 0, {Md - Xdl) I ^2 I = 0. (12) 
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According to our analysis above, all the elements of A/„ and Md can be determined except two of 
them. Suppose that we input the values of these two free parameters, then we can determine 
and Md- Mu and Md can be diagonalized by orthogonal transformation 

V^MuVu = diag(m„, mc, mj), Vj MdVd = diag(mrf, ms,mb). (13) 

By Eq. (fTT|) . Eq. (|13p can be rewritten as 

PuMuPIVu = diag(m„, mc, m*), Vj PdMdP\Vd = diag(md, ms,mb). (14) 
The CKM matrix can be defined as 

VcKM = V^PuP^^Vd. (15) 

By our analysis above, P = PuPd — diag(l, exp'^, exp*''). Therefore, we have four free parame- 
ters, i.e., F„ and Fd respectively in Af„ and M^, and ^ and r\ in P. It is well known that the CKM 
matrix have four free parameters. Hence in principle it is possible that we can adjust the values of 
our free parameters to make them consistent with the experimental data. The similar procedure 
also applies to the lepton sector. 



III. THE APPLICATION TO THE LEPTON SECTOR 

In recent years, neutrino physics has made great progress. The mixing of neutrinos has been 
confirmed and the structure of neutrino mixing matrix has been determined to a reasonable precision 
[T] . It is well known that the neutrino mixing matrix can be expressed with the tribimaximal matrix 
approximately [H]. It is 







V6/3 \/3/3 
-V6/6 \/3/3 V2/2 
\/6/6 -\/3/3 72/2 



(16) 



The important step of our method is to choose the eigenvectors appropriately. In order to choose 
eigenvectors in the lepton sector, we notice some investigations below. 
In the lepton sector, Koide [6] ever suggested an accurate formula 



1/2 V3-\/"^e + m-/^ + m-y 
Foot [7] gave a geometrical interpretation to it, 

(^/m^, Jm^,, ^/m7)(l, 1, 1) 



I (V^' V^) II (1' 1' 1) I ' 



(17) 



(18) 



where (-^ 



/m7) and (1,1,1) are interpreted as two vectors, and 9 is the angle between 



them. If we choose 6* = we get the Koide's mass formula. We have seen that in the tribimaximal 
matrix, the entry 33 is approximate to be Therefore there is a natural connection between 
them. The Koide's mass formula and Foot's geometrical interpretation provide hints for choosing 
the eigenvectors. 

We choose ( y^m^, y^m^, y^ttv) and (1,1,1) as the eigenvectors we want. It is proper that we 
choose (1,1,1) as the eigenvector of neutrino mass matrix belonging to 777,3 (the mass of a neutrino) 
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and (^me, ^ 



/m7) as the eigenvector of the charged lepton matrix belonging to m,-. This 



implies that we designate the vectors as follows 



Vmns = y^jVv 



-k 
-k 



-k 
-k 



P 



\ ^m^+m^+mr ^nic+m^+mr ^m^+m^^+m-r 



2' 
± , 



(19) 



where P = diag(l, exp*^, exp*''). 

As we described above, we need the values of the mass parameters. For the charged leptons, 
the masses are known accurately. However, for the neutrinos, only the mass squared differences are 
measured. The results of global analysis read [U [T^] 



Am?2 = ™2 -m\ = (7.9 ± 0.4) x 10"^ eV^ 
= (2.4 ±0.3) X 10"^ eV^ 



Arn32 1 = 1 TO3 - m\ 



(la) 
(Id) 



If we assume the normal mass hierarchy, i.e., 
then we have 



"^21 




mil 


mis 


7712 




7773 


V "'ii 



0.033 ±0.004, 



= 0.18 ±0.01. 



(20) 
(21) 



(22) 

(23) 
(24) 



It is obvious that the mixing matrix only depends on the mass ratio (TH]. So we only need the 
ratios of neutrinos masses. It is convenient to normalize the neutrinos masses as follows 



Xi = nil = 0.17)7-3, A2 = 7772 = O.27713, A3 = 7)7,3. 



(25) 



Note that these values of mass parameters do not stand for the absolute mass, but just stand for 
the mass ratio. These ratios are consistent with the experimental data we have displayed above. 

Let x = 1, 7/ = 1 and F = 0.279853 ins. in the neutrino mass matrix, by Eqs. ©, ©, 
and (|15p in Appendix C and Appendix D. we can get the neutrino mass matrix 



= 1713 

For the charged leptons, let 

Ai = 77ie = 0.511 MeV, A2 



0.463783 0.279853 0.256364 
0.279853 0.406364 0.313783 
0.256364 0.313783 0.429853 



and 



y = 



105.658 MeV, A3 



:, F = 55.6898 MeV, 



1776.97 MeV, 



(26) 



(27) 



(28) 



similarly we can get the charged leptons mass matrix 

68.2912 55.6898 15.3959 
Ml = I 55.6898 135.509 399.315 
15.3959 399.315 1679.34 



MeV. 



(29) 
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Note that in the course of getting the mass matrices Mi and M^, above, we have chosen the 
eigenvectors of Mi and to be real, hence the mass matrices are real matrices according to our 
analysis in Sec. II. Here we emphasize that this just is a conventional choice, which simplifies the 
equations effectively. 

The matrices that diagonalizes Mi and M^ are given as 

0.599741 -0.800024 0.0164729 \ / 0.169807 -0.798644 0.57735 

Vi ^ I -0.779697 -0.579625 0.23687 , K = -0.776549 0.252265 0.57735 | . (30) 
0.179954 0.154905 0.971402 / \ 0.606743 0.546379 0.57735 

Obviously the third column of Vi is { , , , =, , y. and the third 

^^■m^+m^+rn-r ^m^+m^+m^ ^m^+m^+m^ ' 

column of K is (^, ^, ^)^. 
The MNS matrix is given as 

/ 1 \ / exp»"i \ 
Vmns = V;^PiPIv,P = V;^ \ exp»« \vA exp'"^ . (31) 

\ exp^'' / V 1/ 

The phases ai and 0:2, known as the Majorana phases, have physical consequences only if 
neutrinos are Majorana particles. Because there is no clear evidence whether there is CP-violation 
in the lepton sector, we can not fix the values of ^ and 77. Once we can measure the MNS matrix 
accurately, we can adjust the free parameters F, ^ and 77 to fit the experimental data. If we let 
^ = and 77 = 0, we obtain 

/ 0.816499 -0.577347 \ / exp" 

Vmns = 0.408246 0.577352 -0.707107 exp^"^ | . (32) 

\ 0.408247 0.577352 0.707107 ) \ 

The magnitude of the elements are given as 

/ 0.816499 0.577347 
I Vmns h 0.408246 0.577352 0.707107 | , (33) 
\ 0.408247 0.577352 0.707107 

which is very close to the tribimaximal matrix. 

IV. THE APPLICATION TO THE QUARK SECTOR 

The quark mixing matrix has been determined in high precision. The CKM matrix elements can 
be most precisely determined by a global fit that uses all available measurements and imposes the 
standard model constraints. The allowed ranges of the magnitudes of all CKM elements are [H] 

0.97360 to 0.97407 0.2262 to 0.2282 0.00387 to 0.00405 
VcKM^ vl^Vdh ^ ( 0.2261 to 0.2281 0.97272 to 0.97320 0.04141 to 0.04231 ) . (34) 
0.0075 to 0.00846 0.04083 to 0.04173 0.999096 to 0.999134 

Just like that in the lepton sector, we also need two vectors. We find that a numerical relation 
is well satisfied. It reads 

cos d = ' \, = 0.999549. (35) 

I (^m;;, ^TO^, ^/m7) II [md,ms,mb) \ 

The reasons that we choose this formula are displayed in Appendix A. 
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Hence we might speculate that = 0, while in the Koide's mass relation — j. Because the 
elements of Vckm are given at the scale fi = Mz, we use the quark mass given at the scale fj, — Mz- 
This well satisfied numerical relation suggests us to designate the vectors as follows 



Vc 



KM 



•k 
★ 



P 



y/'niu-\-mc-\-'>Tit \^fnu-\- fn c-i-rrii y/ +mc 



,(36) 



where P — diag(l, exp'^, exp*''). 
It implies that the vector i—r 



is the eigenvector of 

belonging to its eigenvalue nit, and the vector ( , , =, , „ *"% =, , „ =)'^ is the 

eigenvector of Md belonging to its eigenvalue mi,. As we have emphasized before, if we fix the value 
of the free parameter, other elements of the mass matrix are determined. The equations can be 
solved analytically. We display the result in Appendix C. 
For the up-quark sector, let 



Ai = mu{Mz) = 2.33 MeV, As = mc(Mz) = 677 MeV, A3 = mt{Mz) = 181000 MeV, (37) 



and 



/nil, 



V 



/mt \/iTit 
we can get the up-quark mass matrix by Eqs. ([S 



F = 133.18 MeV, 



®, (Uni), (HSl) and (Ull) in Appendix C. 



(38) 



17.9519 133.18 641.198 
Mu = ( 133.18 1335.65 10987.5 
641.198 10987.5 180326 



MeV. 



For the down-quark sector, let 
Ai = nidiMz) = 4.69 MeV, A2 = m,{Mz) = 93.4 MeV, A3 = mb{Mz) = 3000 MeV, 



and 



^^md^ m,^ ^ = 7.99 MeV, 
mt nib 

similarly we can get the down-quark mass matrix 

5.39717 7.99 4.43281 
Afrf = I 7.99 95.5146 90.4138 
4.43281 90.4138 2997.18 



MeV. 



(39) 

(40) 
(41) 

(42) 



Again, in the course of getting the mass matrices M„ and M^, for simplicity, we choose the 
eigenvectors to be real, hence the mass matrices are real matrices by the analysis in Sec. II. 
The matrices that diagonalize M„ and Md are given as 



0.990087 -0.140407 0.00358117 
K = ( -0.140363 -0.988217 0.0610438 
0.00503202 0.0609414 0.998129 



, Vd = 



0.996046 
-0.0888268 
0.00120832 



-0.0888212 0.00156257 
-0.995561 0.0311182 
0.031134 0.999514 



(43) 
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Obviously the third column of Vu is ( , , , , , , . v"^* \t ^j^g third col- 



The CKM matrix is given as 



umn of V, is ^ , , J 

-y/m^+7n|+m^ -v/m^+mj+m^ A/m^+m|+m^ 



/ 1 \ 
VcKM = VjP^P^^Vd = Vj \ exp*« Vd. (44) 

\ exp^'' / 

Given ^ — O.TGtt + j-^ys ^nd rj = O.TGtt, the CKM matrix equals to 

/ 0.973722 + 0.000729014i -0.227558 + 0.0G823131i 0.00224113 + 0.00318905i \ 
VcKM = I -0.227538 + 0.00515369i -0.971078 + 0.0584958i -0.0139226 + 0.0399092i .(45) 
\ 0.00810099 + 0.000510366i 0.0376987+ 0.0177396i -0.729142 + 0.683045i / 

The magnitude of the elements are 

/ 0.973722 0.227707 0.00389778 \ 
I VcKM h 0.227596 0.972838 0.042268 . (46) 
\ 0.00811705 0.041664 0.999099 / 

The quantities of rephasing invariance are calculated as 
a = 02 = arg[-^^^] = 101.59°,/? = = arg[-^^] = 22.74°, 7 = 03 = arg[-|^] = 55.67°, (47) 

VudVy^b ^tdVtb VcdVcb 

and the invariant measure of CP violation is calculated as 

J = -MKdV^fcKTfcKd) = 3.01513 X 10-^ (48) 
while the the experimental data [14j are given as 

a^imtfr, /3= (21.70™)°, 7=(63tl^)°, J = (3.08t°:16) x 10"^ (49) 
They are consistent with each other. 



V. CONCLUSION 

We have illustrated the method in Sec. II, and apply it to the lepton sector in Sec. Ill and to the 
quark sector in Sec. IV respectively. The character of this method is to designate the eigenvector and 
the eigenvalue for the mass matrix appropriately. In the lepton sector, we use the Koide's formula, 
and in the quark sector, we use a similar formula that is well satisfied. Now we give some comments 
about the mass formula we used. 

(1) In the lepton sector, the Koide's mass formula, 

1 _ ^/m^ + ^/fn^ + ^/to7 ^^^^ 
•\/2 \/3 V"^e + Tnn^ + nir 

is satisfied in high precision. It is energy scale insensitive 15 , and its other characters were also 
discussed [ 16j . Several explanations that can realize this formula have been given 17] . Among these 
explanations. Foot's geometrical interpretation seems fascinating phenomcnologically. However, 
there is still no a theoretical model that can realize it. It seems that our method can implement this 
interpretation. If there is no CP-violation and the MNS matrix is tribimaximal, the righthand of the 
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mass formula is connected to the element of the MNS matrix by Foot's geometrical interpretation. 
This is an approach that can lead to the Koide's mass formula. 

(2) In the quark sector, the mass formula like Koide's is unsuccessful [TH]. However, because 
the CKM matrix is very different from the MNS matrix, such a mass formula is useless for us. 
Alternately, we find another well satisfied mass formula. 



This mass formula is connected to the element of the CKM matrix. 

(3) There are two reasons that we choose the mass formula as the vectors. First, the vectors 
are expressed by the mass parameters, so we do not need to introduce extra parameters. Obviously 
that this is an economic choice. Second, there exists such mass formula. They are excellent and 
are well satisfied in high precision, like the Koide's mass formula, but we can not realize them in a 
concise and convincing way. Our method provides an approach that can realize them, nevertheless 
in the special situation if there is no CP-violation. Of course other choices of the vectors are also 
permitted. 

Finally we give some comments about the texture zero structure and our method. The differences 
between our method and the texture zero structure are: in our present approach, at first we choose 
one eigenvector of the mass matrix, and then we can determine other elements of the mass matrix; in 
the texture zero structure, some elements of the mass matrix are supposed to zero, which equals to 
designate the eigenvectors, but we do not know the eigenvector at first. Therefore, in our approach, 
we have the freedom to choose the eigenvectors to satisfy other request. For example, the Koide's 
mass formula can be realized in our approach through Foot's geometrical interpretation. 

As the texture zero structure, our approach is also consistent with current experimental data, 
and this approach has the merit that it can realize some well satisfied mass formula, for example, 
the Koide's mass formula and the mass formula suggested by us. However, just as many texture 
zero structures, our approach is purely phenomenological, and there is still no a theoretical model 
to realize it. Therefore it is worthy to investigate whether our approach can be realized in some 
theoretical models. If this is true, it will provide a new theoretical and phenomenological approach 
to deal with the mass and mixing problems. Besides, we point out that in our paper we just restrict 
our discussions in a simple case, in which we choose the eigenvectors to be real, hence we get the 
real mass matrices. This is just a conventional choice that simplifies the equations. Other cases, 
in which the eigenvectors are complex, are permitted, and they should be considered if they are 
needed by some underlying theories. 
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APPENDIX A 

Table: Quark masses at the Z mass scale in the standard model [H] 




= 0.999549. 



(51) 




-0.11 



,+U.tiU 



McV 
MeV 
GeV 



With the values of the quark masses above. 



TO„(Mz) = 2.33 MeV, mdMz) = 677 MeV, mt{Mz) = 181 GeV, 
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md{Mz) = 4.69 MeV, ms{Mz) = 93.4 MeV, mb{Mz) = 3.00 GeV, 



cos6' 



( ^/rnl) {md,ms,mb) 

(v^'V^'V^) II {md,ms,mb) 

cos 6* = 0.999549 is very close to 1. 

Note that we do not choose the mass formula below, 



0.999549, 



cos^ 



COS 9 



cos 9 



I {mu,mc,mt) \ \ y/Ws, y/mZ) \ 

_ [mu ,mc,mt) {md ,ms,mb) _ 
I {mu,mc,mt) \\ {md,ms,mb) \ 



= 0.984685, 

0.999624, 

0.992936. 



As we have shown in Sec.IV, the mass formula is connected to the element of the CKM matrix by, 

fnu\/Tnd + exp'^ mc^Prnl + exp*'' mt^Jrub 



V33 



V33 



mumd + exp*f rricms + exp"' mtrrib 



V3s = 



i/m2 +ml+ m^^m'j + ml+ml 

'm7i^~ + exi)'^ ^7777:^7777 + oxp"' ^7777 



■v/to-u + mc + mt\/md + ms+mb 



and 



rriuJmd + exp*^ rricJrns + exp*** mtJmb 



0.983385 ^1 



0.999391 ^ 



"v/m^ + + mly/rrid + ms+mb 
murud + exp'^ rricms + exp"* mtrrib 



1^ 0.984685, 



^ 0.999624, 



rriu\/md + exp'^ JfricJfris + exp"' JmtJrrib 



0.971465 <| 



1^ 0.992936, 



i/m„ + mc + mt\/md + rog + mft 
but the experimental value of | V33 | is 

0.999096 <| V33 |< 0.999134. 

They are not consistent. Therefore, in order to make them consistent with the experimental data, 

wc must change the mass parameters that wc have \iscd. But the mass formula we have chosen does 
not have this problem. Our choice is just a convenient one. Of course other choices are permitted 
if they are consistent with the experimental data. 



APPENDIX B 

As we have emphasized in the text, the analytical expressions of the solutions will be complicated. 
In this appendix, we give some analysis that will simplify the expressions effectively. 
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In order to express the other elements of the mass matrix with the quarks masses and the free 
parameter F, we have to solve the equations displayed below 







( A \Wi -U P^r-ri -U D'To — 

l^/l — A)Jb\ -T X2 — 


n 






Pt, -4- 


( U \ Wo -I- P nvn^" T-o — 


n 


('9'! 


iA - \A{B - \AiC - \^) - E'^iA- 

^ ^*i7\-'--' -1-/ V 


-All - D'^iB - \^)- 


F'^{C - Ai) + 2DEFN = 







[A - \2){B - M){C - A2) - E\A - 


- A2)-Z?2(B-A2)- 


F^{C - \2) + 2DEFN = 


0, 


(4) 


{A-\:i){B-\^){C-\3)-E\A 


-A3)--D'(B- A3)- 


- F^{C - X3) + 2DEFN ^ 


: 0. 


(5) 


It is well known that 










Trace(M) = Ai 


+ A2 + A3 =^ A + S 


+ C = s, 




(6) 


in which s = Ai + A2 + A3. 










If Ai ^ A2, by Eqs. @, Q and (HD, 


we obtain 








AB + BC + 


AC + {D^ ^E"^ - 






(7) 



in which r = A1A3 + A1A2 + A2A3. 

By Eqs. ([T]) and ([2]), we can express D and E in terms of A, B and F 

£» = -{A-\)x - Fy, 
E = [-Fx -{B- \)y] exp~'", 

in which we have let a; = ^ and y — supposing that 2:3 ^ 0. 
We obtain five new equations 



{A - A3)(B - A3)(C - A3) - E^A - A3) -D\B- A3) - F\C - A3) + 2DEFN ^ 0, 

A + B + C ^s, 
+ {D^ + E^ + F^), 
D = -{A - X)x~ Fy, 
E = [-Fx -{B- \)y] exp-^" . 



AB + BC + AC 



(8) 
(9) 



(5) 
(6) 
(7) 
(8) 
(9) 



Let N = cos a = ±1, by these equations, we can express A, B, C, D and E in terms of F. 

In the following we will deduce another formula, which will simplify the expressions effectively. 
By Eq. ([6|), we obtain 



C = s-A-B. 

We have argued that N = cos a = ±1 in Sec. II. So by Eq. ©, we have 

E ^ ±[-Fx- {B - \)y]. 

In Eq. ([7]), we express C, D and E in terms of A, B and F by Eqs. (H 
some simplification, we obtain 



aB^ + bB 



0, 



(10) 
(11) 

(Uni) and (HI]). After 
(12) 



in which 



a = 1 + r, 

b = A- s- 2Xy^ + 2Fxy, 

c={\y- Fxf + [{A - \)x + Fyf + F^ + r - A{s - A). 
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Then B can be solved as 



^^-b±VlZ^^ (13) 
Hence if we know the values of A and F, the value of B is determined by the Eq. ([T^ . 
APPENDIX C 

In this appendix, we give the analytical expressions of the solutions for the quark sector. 

In the quark sector, let A = A3 and a = 0. According to Eqs. ([U, ([S]) and dH), we can express 
_B, D and E with A, C and F linearly. The Eqs. ([5|) and ([7]) are left invariant. We express B, D 
and E in terms oi A, C and F in Eqs. ^ and d?]). Then in Eqs. ^ and jT]), only ^, C and F 
are present. Therefore by Eqs. ([5]) and ([7]), we can solve A and C in terms of F. The solution is 
displayed as follows 



in which 



a' = x"* + x^y^ + 2j:^ + + 1, 

6' = -2A3J;'' + 2Fyx^ - (Ai + A2 + 2A3y2 + 2A3)x2 + {2Fy^ + 2Fy)x- 
(l + y2)(Ai+A2), 

c' = A^x^ - 2A3^^2/x3 + + + A^yS + A1A3 + A2A3)a;2 - ^X^Fiy^ + y)x^ 
F\y^ + 2y^ + l) + {y^ + l)\^\2, 

The other elements of the mass matrix are given as 



- - ^i^^, (13) 

C = s-A-B, (10) 

D = -iA-X3)x~Fy, (8) 

E = -Fx-iB-X3)y. (9) 



APPENDIX D 

In this appendix, we give the analytical expressions of the solutions for the lepton sector. 
In the lepton sector, the analysis in Appendix B and C applies similarly. Let A = A3 and a — 0. 
We find that the solutions we need can be analytically expressed as follows 



A^±±^^I^, (15) 

in which b' , a', c' are expressed as the same as that in Appendix C. Eqs. ([S]), (O, pU|) and in 
Appendix C apply similarly. 
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